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We have modeled the time-variable profiles of the Ha emission from the nonaxisymmetric disk and debris tail 
created in the tidal disruption of a solar-type star by a 10 6 Mq black hole. Two tidal disruption events were 
simulated using a three dimensional relativistic smoothed particle hydrodynamic code to describe the early 
evolution of the debris during the first 50-90 days. We have calculated the physical conditions and radiative 
processes in the debris using the photoionization code CLOUDY. We model the emission-line profiles in the 
period immediately after the accretion rate onto the black hole becomes significant. We find that the line profiles 
at these very early stages of the evolution of the postdisruption debris do not resemble the double-peaked profiles 
expected from a rotating disk, since the debris has not yet settled into such a stable structure. As a result of 
the uneven distribution of the debris and the existence of a "tidal tail" (the stream of returning debris), the line 
profiles depend sensitively on the orientation of the tail relative to the line of sight. Moreover, the predicted 
line profiles vary on fairly short time scales (of order hours to days). Given the accretion rate onto the black 
hole we also model the Ha light curve from the debris. 



1. INTRODUCTION 



A star in an orbit around a massive black hole can 
get tidally disrupted during its close passage by the 
black hole. After several orbital periods the debris 
from the disrupted star settles into an accretion disk 
and gradually falls into the black hole [E H H 13 As 
material gets swallowed by the black hole intense UV 
or soft-X ray radiation is expected to emerge from the 
innermost rings of the accretion disk H HH IE H EH 
[H E2. For black hole masses M bh < 10 7 Mq, tidal 
disruption theory predicts flares with luminosities of 
the order of the Eddington luminosity with durations 
of the order of months, and with spectra that peak in 
the UV/X-ray band [J El El El El High-energy 
flares from the central source illuminate the debris, 
the photons get absorbed, and some are re-emitted 
in the optical part of the spectrum (i.e. the light is 
"reprocessed" ). One of the spectral lines in which this 
phenomenon can be observed is the Balmer series Ha 
line (X r est = 6563 A). 

The disruption of a star begins when the star ap- 
proaches the tidal radius, r t ~ r^Mbh/M*) 1 ' 3 , the 
point where the surface gravity of the star equals the 
tidal acceleration from the black hole across the diam- 
eter of the star (r* and M* are the radius and mass 
of the star and M^h is the mass of the black hole). A 
10 6 M Q black hole is often used as a prototypical ex- 
ample in tidal disruption calculations. This choice is 
motivated by the criterion for a solar-type star to be 
disrupted before it crosses the black hole event horizon 
(i.e. the Schwarzchild radius, r s ) in order for emission 
to be observable. For supermassive black holes with 
Mi,h > 10 8 Mq, r s > r t and the star falls into the 
black hole before it gets disrupted. 



2. CALCULATION 



Tidal disruption simulations were carried out with 
a three dimensional relativistic smoothed particle hy- 
drodynamics (hereafter SPH) code in order to study 
the dynamical evolution of the post-disruption debris. 
The SPH code used provides a description of relativis- 
tic fluid flows in a static curved spacetime geometry 
[HI Il8| . We use it to simulate the tidal disruption of 
a star in the potential of a Schwarzchild black hole. 

Two different simulations were carried out with 
5,000 and 20,000 particles (hereafter 5k and 20k re- 
spectively) contributing equally to the mass of a 1 M Q 
star. The 5k calculation follows the debris for 94 days 
in total. After 34 days, significant accretion onto the 
black hole begins. Our investigation follows the evo- 
lution of the line profiles in the last 60 days. The 20k 
simulation spans 53 days, during which the evolution 
of the line profiles was followed for the last 6 days. Us- 
ing both the 5k and 20k simulations in the line profile 
modeling we take advantage of the longer time span 
in the former and better resolution achieved with the 
larger number of particles in the latter. Figure^shows 
particle distribution maps after the second pericentric 
passage, at the beginning of the accretion phase and 
at the end of the 5k simulation. At the early stages 
of the tidal event most of the particles were located in 
the pronounced tidal tail. Sixty days later, about 20% 
of the particles are scattered from the tidal tail and 
form a quasi-spherical distribution, with most of its 
mass concentrated in the equatorial plane. We refer 
to the spheroidal part of the debris as the halo and to 
its planar component as the disk. There is a concern 
that some fraction of particles of the halo are a nu- 
merical artifact of the simulation. We further mention 
the implications of existance of the spheroidal halo for 
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the total Ha luminosity and emission line profiles in 

inland mi 

We follow the line profile calculations carried out 
by ^3 and j2£j to obtain the observed profile from 
a Keplerian, relativistic, thin disk in the weak field 
approximation. The main objective of the calculation 
is to obtain the final expression for the flux density 
in the observer's frame as a function of parameters 
defined in the reference frame of the debris. 

In order to calculate the resulting Ha luminosity of 
the debris it is necessary to determine the efficiency 
with which the debris reprocesses the incident radia- 
tion. This efficiency can be characterized by the sur- 
face emissivity of the debris as a function of radius, 
e = €q where eo is a constant and £ = r/r g is the 
dimensionless radius (r g — r s /2 = GM^jc? = Mi,/,,). 
We have used the photoionization code CLOUDY [2l| 
for numerical calculations of the response of the de- 
bris to illumination. For a more detailed description 
of the calculation refer to l22l. 



3. RESULTS 
3.1. Light Curves 

With the above emissivity prescriptions we have cal- 
culated the observed Ha luminosity curve of the de- 
bris at a particular time step by computing the time 
at which the light was emitted from the debris and 
by finding the ionizing flux that was illuminating that 
location at the time the emission occurred, according 
to the light-travel time from the black hole to that 
particular region of the debris. We assumed that the 
observer is located in the x'z' plane, above the orbital 
plane at i = 30° to the z'-axis, at a distance d — > +oo. 

Figure [21 shows three different Ha light curves from 
the debris confined to a plane (assuming £j„ = 500, 
£,out = 40, 000) during the 60-day accretion phase of 
the 5k simulation. Figure^ shows the accretion lumi- 
nosity on a logarithmic scale (solid curve), calculated 
from the accretion rate of the debris in the simula- 
tion. The UV/X-ray luminosity curve is arbitrarily 
scaled and overplotted on the top of the Ha curve for 
comparison. The Ha light curve departs from the ac- 
cretion light curve at late times, though the departure 
appears small in the logarithmic plot (used here due 
to the large dynamic range of the light curves). The 
same effect is more noticeable in Figure 03, where the 
accretion luminosity is proportional to i~ 5 / 3 and the 
Ha light curve is plotted on a linear scale. The Ha 
light curve roughly follows the shape of the incident 
UV/X-ray light curve at early times but decays faster 
at late times. The faster decay in the Ha light curve 
reflects the debris evolution in time: as the tail be- 
comes more elongated, the incident photons travel a 
longer way to illuminate the debris. Consequently, 
the intensity of the illuminating light gets lower in 
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Figure 1: Maps showing the positions of the SPH 
particles from the 5k simulation after the second 
pericentric passage, projected in the x'y'-plane at two 
different times. Particles gravitationally bound to the 
black hole are colored black, while unbound particles are 
colored grey. The dashed line represents the initial 
trajectory of the star before disruption and the 
trajectory of the center of mass of the debris after 
disruption. The maximum particle velocities are of order 
10 _2 c. Upper Panel: Particle map at the start of the 
accretion phase, 34 days after the disruption occurred. 
Inset: Particles in the inner region of the debris, orbiting 
close to the black hole. Lower Panel: Particle map at the 
end of the simulation, 94 days after the disruption. The 
tidal tail can be clearly separated into particles that are 
unbound and about to escape the black hole and 
particles that are returning towards the black hole. The 
inner region of the debris consists of returning particles 
from the inner tail that have been scattered and form a 
flat rotating structure around the black hole. 



the later stages of the tidal disruption event. At late 
times the decay of Ha light curve stops because of the 
return of more particles from the tail to the imme- 
diate vicinity of the black hole. As particles diffuse 
from the high density tail to the lower density disk, in 
later stages of the simulation, their emission efficiency 
increases and they contribute a significant amount of 
Ha light to the light curve. To isolate the effect of the 
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debris evolution in time from the evolution of the illu- 
minating light curve, we calculate the Ha light curve 
in the case of constant illumination (Fig. 120 ■ Here, 
the relative departure of the Ha light curve from the 
UV/X-ray light curve can be interpreted as a con- 
sequence of the expansion and redistribution of the 
debris. The Ha luminosity appears to level off at late 
times because the debris disk begins to settle into a 
quasi-steady configuration. 
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Figure 2: The Ha light curves (triangles) resulting from 
the reprocessing of three different illumination light 
curves. The solid line shows a scaled version of the 
UV/X-ray continuum light curve that illuminates the 
debris, which in (a) follows the SPH accretion rate, (b) 
decays with time as t~ 5 ^ 3 , and (c) is constant. Note that 
the vertical axis is calibrated logarithmically in (a), while 
in (b) and (c) they have the same linear scale. The Ha 
light curves rise initially as the illumination front 
propagates through the debris and then decay faster than 
the UV/X-ray light curves. 

The observed Ha flux depends sensitively on the 
UV/X-ray light curve, on the distribution of matter 
that makes up the inner portion of the debris, and on 
how quickly particles redistribute themselves in phase 
space. The main features of the Ha light curve are: 
an initial rise followed by a decline, with superposed 
fluctuations. The initial rise is a consequence of the 
propagation of the initial illumination front through 
the debris. The fluctuations are a result of the fluc- 



tuations in the accretion rate, which are caused, in 
turn, by the finite number of particles employed in 
the simulation. The decay rate of the Ha light curve 
is determined by the decay rate of the UV/X-ray light 
curve, debris expansion and redistribution rate. 

The CLOUDY calculation predicts a time-average 
Ha luminosity from the tidal debris of about 10 36 , 
10 37 , and 6 x 10 38 erg s _1 for the tail, disk and 
halo, respectively. In the earlier stages of the disrup- 
tion event when the UV/X-ray luminosity is super- 
Eddington, the Ha luminosity is expected to be up 
to 80 times higher than its average value and compa- 
rable to that of tidal disruption candidates observed 
in the local universe (see Figure [21 ■ These examples 
include NGC 4450 (at 16.8 Mpc) with an Ha lumi- 
nosity of L Ha =1.8xl0 39 erg s" 1 and NGC 1097 
(at 22 Mpc) with L ffQ =7.7xl0 39 erg s" 1 |3. Thus 
the emission-line signature of a tidal disruption event 
should be detectable at least out to the distance of the 
Virgo cluster. In practice, however, the detection of 
such emission lines from low luminosity sources may 
be complicated by their weak contrast relative to the 
underlying stellar continuum. 

3.2. Line Profiles from the Debris 

In Figures EHHl we show sample line profiles to il- 
lustrate how they evolve in time and how they are 
affected by the choice of model parameters and by 
the orientation of the observer. Figure is a "trailed 
spectrogram" summarizing the temporal evolution of 
the line profiles from the two different SPH runs; it 
is a 2-dimensional map of the Ha emission as a func- 
tion of projected velocity and time. Figure 0] shows 
a different representation of the evolution of the line 
profile with time, which effectively comprises selected 
time slices from the trailed spectrogram. Figures |S] 
and El show how the inner radius of the line-emitting 
region and the azimuthal orientation of the observer 
affect the observed line profiles. 

A property that is immediately obvious in the line 
sequence is the change of the profile shape with time 

(Figures and . The adopted low value of ve- 
locity dispersion allows us to resolve individual parti- 
cles in the trailed spectrograms, orbiting around the 
black hole. The evolution of the line intensities in 
time roughly follows the behavior of the UV/X-ray 
luminosity but decays somewhat faster in time. The 
multi-peaked line profile is a consequence of the ve- 
locity field of the inner debris, which consists of the 
inner portion of the tidal tail that is falling towards 
the black hole (towards the observer) and debris that 
is rotating around the black hole after being scattered. 
The line profiles and their variability could be obser- 
vationally important features of the debris just formed 
from tidal disruption. The variable line profiles might 
be observed and recognized on the relatively short 
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Figure 3: Trailed spectrogram of the simulated Ha 
emission-line profiles from the 5k simulation spanning 60 
days (left) and from the 20k simulation spanning 6 days 
(right). This is effectively a 2-dimensional intensity map 
versus projected velocity of the emitting material and 
time. Darker shades correspond to higher intensity. The 
scale on the right represents time since the tidal 
disruption event. 



time scale of hours to days. 

The profiles become broader as the inner radius 
of the line-emitting regions decreases since higher- 
velocity gas resides at smaller radii (see Figure EJ. 
The approximate full width at zero intensity of the 
profiles ranges from 4, 500 km s _1 for = 10, 000 to 
18,000 km s _1 for = 200. We find that line pro- 
files change from the profiles dominated by the emis- 
sion red-ward from the rest wavelength for £j„ < 1000 
to narrower profiles dominated by the blue- ward emis- 
sion from the tail for £j„ > 1000, since for large values 
of £m, the high- velocity rotating gas in the vicinity 
of the black hole is excluded and the dominant con- 
tributions to the line profile come from the tidal tail. 
The intensity of the line also decreases with increasing 
inner radius, making the outer regions of the debris 
harder to observe. 

Because of the non-axisymmetric geometry and ve- 
locity field, the line profiles emitted by the debris, de- 
pend on the orientation of the tidal tail relative to the 
observer. In Figure El we show the effect of azimuthal 
orientation (f>o of the debris, with respect to the ob- 
server. The values of O are 45°, 90°, 120°, 180°, 
220° and 270°, as measured in a counterclockwise di- 
rection from positive x'-axis to the observer's line of 
sight. These can be compared with the profile corre- 
sponding to the same time in Figure for <fio = 0° . 
The position of the peaks in Figure varies relative 
to the rest wavelength, since the relative direction of 
bulk motion of the material depends on the observer's 
orientation. For example, it is possible to distinguish 





=500-40000 
=30" 










r=100 km a 


W 




6 d g h oo m 




3 d 18 h 00™ 






6 d 6 h 00 m 




3 d 10" 30 m 






8 d 2 h 15 m 




3 d 3 h 00™ 


m 




5 d 22 h 30 m 


.^M 


2 d 19 h 30 ra 






5 d 18 h 45 m 


■Jli+t, Ah 12* °°~ 






5 d 15 h 00™ 


.,klV\n SUN z a 4" 30- 






5 d ll h 15 m " 










5 d 7 h 30 m 


a A%J>w"\JL 1* I* a °" 






5 d 3 h 45 m 




jd 6 h 00 m 






5 d 00 h 00 m 




d 22 h 30 m " 


in 
d 


, jwAjlJ 


4 d 20 h 15 m ~ 


.A*. A 


d 15 h 00 ra 




, ^aK Jl 


4 d 16 h 30 m 




d 7 h 30™ 




...aA /u. 


4 d 13 h 45 m 


, Aj, .,. nfl 


Relative time 
d l h 48™ 






Relative time 
4 d 9 h 00 m 



6200 6400 6600 6800 7000 6200 6400 6600 6800 7000 
Wavelength (A) Wavelength (A) 



Figure 4: Sequence of Ha profiles emitted from the 
region £ £ (500, 40 000) over a period of 6 days (20k run). 
The relative time from the beginning of the accretion 
phase onto the black hole is marked next to each profile. 
The accretion phase begins 47 days after the tidal 
disruption. The inclination of the debris plane and the 
velocity shear are as marked on the figure. 

the emission from the tail for the range of azimuthal 
orientations 90° — 220°. The tail emission in these 
profiles appears as the most blueshifted peak, since 
these are the orientations for which different portions 
of the tail flow towards the observer. 

We have computed model profiles for several dif- 
ferent X-ray illumination light curves keeping all the 
other parameters fixed. We used (a) the light curve 
obtained from the accretion rate in the 5k simulation, 
(b) the light curve from the accretion rate as predicted 
by theory, i.e. oc i" 5/3 0,13, and (c) a light curve 
that is constant in time (Figure [5J. We find that the 
line profile shapes do not depend sensitively on the 
shape of the light curve. This is a consequence of the 
centrally "weighted" emissivity profile of the debris 
which causes the innermost emission region to be the 
dominant contributor of the Ha light. In the inner- 
most region of the debris the dynamic range in light- 
travel times is not large; therefore the illumination of 
the innermost emitting region is almost instantaneous. 
Over the very short light-crossing time of the central 
emitting region, the gradient in the UV/X-ray light 
curve is small and the illumination is nearly constant 
over this region. The fast fluctuations in the illumi- 
nating light curve on the other hand are smoothed 
out during reprocessing in the debris, and cannot be 
identified in the Ha light curve. 

The temporal variability of the Ha emission line 
profiles from the post-disruption debris is one of the 
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Figure 5: Ha emission line profiles simulated for seven 
different values of inner radius (£in), as determined by 
the propagation of the ionization front through the 
debris. The relative time for profile frames is 6 d 6 m . 
The intensity of the profile calculated for £j n = 10, 000 is 
multiplied by the factor of 100. The inclination and 
velocity shear are as marked on the top of the figure. 



Figure 6: Ha emission line profiles simulated for six 
different azimuthal orientations of the debris with respect 
to the observer, as represented by (f>o- See profile in 
Figure |H with the time label 5 d 22 h 30 m for orientation 
(f>o = 0°. The size of the emitting region, inclination and 
velocity shear are as marked on the top of the figure. 



4. CONCLUSIONS 



important indicators of a tidal disruption event. In 
order to capture the rapid profile variability, due to 
the variable illumination, the exposure time should 
be comparable to the light crossing time of the in- 
nermost regions of the line-emitting debris, which has 
the fastest and strongest response to the ionizing ra- 
diation. Longer exposures are expected to capture 
the average shape of the rapidly varying line profiles. 
The light-crossing time of the innermost regions of the 
debris is about 8 £2 M 6 minutes (where £ 2 = £/100 
and M 6 = M^/IO 6 M©), while the exposure times 
are typically about 30-60 minutes (for galaxies at the 
distance of the Virgo cluster, for example). Thus, if 
an event is caught early in its evolution and the light- 
crossing time is relatively long (i.e., M > 10 6 M Q ), 
there is a chance of detecting variability caused by 
changing illumination over the course of one to a few 
nights. On longer time scales, variability is caused by 
changes in the structure of the debris. In the presence 
of the spheroidal halo, the variability of the lines may 
be modified by the long diffusion time scale of photons 
through the halo. The component of the tidal tail out- 
side the halo will then still respond to the variability 
but on the time scale set by the light reprocessed by 
the halo. 



We modeled the emission-line luminosity and pro- 
file from the debris released by the tidal disruption 
of a star by a black hole in the early phase of evo- 
lution. Our model predicts prompt optical evolution 
of the post-disruption debris and profile shapes dif- 
ferent from circular and elliptical disk model profiles. 
Since line profiles observed so far in LINERs (low- 
ionization nuclear emission regions) look more disk- 
like and evolve slowly, the observations are likely to 
have caught the event at late times (> 6 months after 
the initial disruption) , after the debris has settled into 
a quasi-stable configuration. 

The line profiles can take a variety of shapes for 
different orientations of the debris tail relative to the 
observer. Due to the very diverse morphology of the 
debris, it is almost impossible to uniquely match the 
multi-peaked profile with the exact emission geome- 
try. Nevertheless, the profile widths and shifts are 
strongly indicative of the velocity distribution and the 
location of matter emitting the bulk of the Her light. 
Profile shapes do not depend sensitively on the shape 
of the light curve of the X-rays illuminating the debris. 
They strongly depend on the distance of the emitting 
material from the central ionizing source, which is a 
consequence of the finite propagation time of the ion- 
ization front and the redistribution of the debris in 
phase space. It may be possible to distinguish be- 
tween the two effects observationally, based on their 
different characteristic time scales. 
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If X-ray flares and the predicted variable profiles 
could be observed from the same object they could be 
used to identify the tidal disruption event in its early 
phase. The X-ray flares can be promptly detected by 
all-sky synoptic X-ray surveys and high energy burst 
alert missions such as Swift. The evolution of the tidal 
event may then be followed with optical telescopes 
from the ground on longer time scales and give an 
insight in the next stage of development of the debris. 
Thus, simulations of the tidal disruption process on 
longer time scales (of order several months to a few 
years) are sorely needed. Calculations of the long- 
term evolution of a tidal disruption event can predict 
the type of structure that the debris finally settles 
into and whether its emission-line signature resembles 
the transient double-peaked lines observed in LINERs. 
This study would provide an important insight into 
the evolution of LINERs. 

Finally, the observed rate of tidally disrupted solar 
type stars can constrain the rate of captured com- 
pact ob ject s which are important gravitational wave 
sources [25|, ul& > an d the capture rate of main se- 
quence stars in our Galaxy, which are expected to emit 
the peak of the gravitational radiation in the LISA fre- 
quency band and can be detected in the local universe 
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